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Abstract 

The quantum resonances of an atom possessing a single valence electron which shows hyper- 
fine interaction with the nucleus is investigated in the presence of a three dimensional magnetic 
quadrupole field. Particular emphasis is put on the study of the interplay of the hyperfine and 
quadrupole forces. Analyzing the underlying Hamiltonian a variety of symmetries are revealed 
which give rise to a two-fold degeneracy of the resonance energies. Our numerical approach em- 
ploys the complex scaling method and a Sturmian basis set. Several regimes and classes of short- 
lived and long-lived resonances are identified. The energies and decay widths of the resonances are 
characterized by their electronic and nuclear spin properties. 
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I. INTRODUCTION 



Ultracold atomic gases offer a wealth of opportunities for studying quantum phenomena 
at mesoscopic and macroscopic scales [l[ 2I, y, |j, 0]. The majority of experiments in the 
ultra-low temperature regime have so far been performed with alkali atoms. Their ground- 
state electronic structure is characterized by the fact that all electrons but one occupy closed 
shells. The latter, the so-called valence electron, is situated in the s-orbital of the outermost 
shell. The coupling of the nuclear spin I to the total angular momentum of the active 
electron via the hyperfine interaction causes a splitting of the electronic energy levels into 
several branches which can be characterized by the quantum number of the total spin F. 
Due to a vanishing orbital angular momentum this hyperfine interaction consist exclusively 
of the coupling of the electronic spin S = |, to the nuclear spin I. This leads to the two 
possibilities F = I ± S. 

Inhomogeneous magnetic field configurations represent a key component for the control 
of the motion of cold atoms, specifically for their trapping. The underlying Zeeman in- 
teraction occurs due to coupling of the total magnetic moment of the atom to the external 
magnetic field. Quantum states and particularly resonances of atoms in inhomogeneous field 
configurations are therefore of immediate interest to cold atomic physics in general and .have 
been investigated in the past by several groups. This includes the quadrupole field 



the wire trap jg, |9|, llO, LL3J and the magnetic guide and Ioffe trap [12|, [13j, UM Il6j |. In 



these traps we typically encounter quantum resonances with a certain lifetime instead of 
stationary 'stable' quantum states. Very recently [it| a new class of short-lived resonances 
possessing negative energies for the case of the 3d magnetic quadrupole field was found. 
These resonances originate from a fundamental symmetry of the underlying Hamiltonian. 
In contrast to the positive energy resonances, they are characterized by the fact that the 
atomic magnetic moment is aligned antiparallel to the local direction of the magnetic field, 
and their lifetimes decrease with increasing total magnetic quantum number. 

In all of the above investigations it is assumed that the hyperfine interaction is much 
stronger than the interaction with the external inhomogeneous magnetic field and therefore 
the electronic and nuclear angular momenta firstly provide a total angular momentum which 
then interacts with the magnetic field. This physical picture truly holds for (alkali) atoms 
in their electronic ground state and macroscopic as well as microscopic (atom chip) gradient 
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fields. The atom is consequently treated as a point particle with the total angular momentum 
F. In the present work we study the case for which both interactions, the hyperfine and the 
field interaction, have to be taken into account on equal level for the description of the neutral 
atoms in the field. This case is primarily of principal interest but is expected to describe the 
magnetized hyperfine properties of systems with very small hyperfine interactions and/or 
strong gradient fields, such as electronically excited atoms. In the latter case an admixture 
of different hyperfine states F due to the field interaction has to be expected. In ref. [3] 
the possibility to tune the hyperfine splitting by dressing the electronic energy levels by 
a microwave has been demonstrated. Hence, by this method it is thinkable to achieve a 
scenario in which the magnetic and the hyperfine-interaction become comparable even for 
ground state atoms. 

We focus on atoms possessing a single active valence electron with spin S — \ and a 
nucleus with spin / = | in a three dimensional magnetic quadrupole field. In detail we 
proceed as follows. Section [Til contains the derivation of the underlying Hamiltonian. Our 
computational method to calculate the resonance energies and lifetime is outlined in section 
IIHI Section [TV] contains a presentation and discussion of our results. We conclude with a 
summary in section IVl 



II. HAMILTONIAN 



Taking into account the hyperfine interaction, the Hamiltonian describing the motion of 
an atom with mass M and electronic and nuclear spin S and J, respectively, reads in a 
magnetic quadrupole field 

H = m + H " F + HB (1) 

where Hhf describes the hyperfine interaction between the outermost single valence electron 
and the nucleus 

H HF = AI ■ S. (2) 



Here A is a constant which for an s-electron is given by 19] 



A = jj0# e #„/i£/i7v|?MO)| 2 (3) 

g e and g n are the g-factors of the electron and the nucleus, respectively, and V's(O) is the 
value of the valence s-electron wave function at the nucleus. u,b = an d I^n — 



are the Bohr magneton for the electron and the proton, respectively. H B accounts for the 
interaction of the magnetic moment of the electron and the nucleus with the magnetic field 



H B = fx B g e (S + al) ■ B (4) 

with 

a = « . (5) 

The vector of the three-dimensional quadrupole magnetic field is given by B(r) = 
b(x, y, — 2z). By substituting equations (2-5) into equation (1) and performing the scale 
transformation 

pi = {2Mbfi B g e )~ 1/3 Pi and x t = {2Mbfi B g e ) 1/3 Xi 

we obtain 

H = ^[P 2 + x(S x + ocI x ) + y{S y + aly) 

-2z(S z + aI z )+[3I-S] (6) 

with 3 = , . where the bar have been omitted. The energy is now measured in units 
of jj(2Mb[i B g e ) 2 / 3 . For convenience and in anticipation of the forthcoming discussion we 
transform the Hamiltonian to a spherical coordinate system, i.e. (x, y, z) (r, 6, <fi). Writing 
the momentum operator explicitly and using atomic units we obtain 

rr 1 r <9 2 2d L 2 

+ r sin 9 sin <p(S y + al y ) — 2r cos 9(S Z + al z ) 

+ /3I-S]. (7) 

Exploring the symmetries of the Hamiltonian we find 16 discrete symmetry operations 
and a continuous symmetry generated by J z , which is the z-component of the total angular 
momentum of the atom. 

The Hamiltonian ([H]) basically differs from the Hamiltonian (4) in ref . [3] by the additional 
hyperfine interaction term. This term does not introduce new symmetries and degeneracies 
except that all symmetry operations found in ref.pj have now to be generalized to the case 
of the additional presence of the nucleus angular momentum. For a discussion of degeneracy 
and symmetries in more details, we refer the reader to ref. 
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III. NUMERICAL APPROACH 



The Hamiltonian ([6]) does not support bound states. Its continuous spectrum is charac- 
terized by resonances which are localized in space (at t — 0). The time evolution of these 
states is given by 

Mt) = e^Vji(O) (8) 

where E is complex 

E = s-i\. (9) 

Here e and T are the energy and decay width of the resonance, respectively. Because of the 
imaginary part — i~ the resonances decay exponentially with a lifetime of r = In order 
to calculate the energies and decay widths of the scattering wave functions we employ the 



complex scaling method (see ref. 20j and references therein). The complex scaling transfor- 
mation only affects the continuum, while the bound states do not change their positions. 
The continuum states are rotated around the corresponding threshold into the lower half 
complex energy plane. Eigenvalues belonging to resonances, once revealed, maintain their 
positions. 

Since resonance states of the complex scaled Hamiltonian are square integrable the linear 
variational principle can be applied. Expanding the wavefunction in a set of basis functions 
and calculating the expansion coefficients results in a large-scale algebraic generalized eigen- 
value equation which can be solved by employing a Krylov space method [2lj]. The basis 
set must be chosen in such a way that the exact wave function can be approximated to a 
sufficient degree of accuracy by as small as possible number of functions. Therefore the form 
of the basis set must be adapted to the geometry and the symmetries of the system. We 
found a so-called Sturmian basis set of the form 

| n,l, mi ,m s ) m = Rf^Yr^-^iO^mrnWirris) (10) 

suitable for our purposes. Here the functions Y" z m are the spherical harmonics. For fixed m 
the linear variational combination of the basis functions ip™ t m . ms , yields, per construction, 
eigenstates of the Hamiltonian and J z simultaneously. For expanding the radial part we take 
the non-orthogonal set of functions 

R«\r) = e-%L n ((r) (11) 
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where L n (x) are the Laguerre polynomials. The free parameter, (, has the dimension of an 
inverse length and can be tuned to improve the convergence behavior in different regions of 
the energy spectrum. It should be chosen such that l/( corresponds to the typical length 
scale of the states to be approximated. Using the basis set (fTUj) all matrix elements of the 
Hamiltonian([7|) can be calculated analytically. 



IV. RESULTS 



Let us now discuss the results we obtained while studying an atom with hyperfine in- 
teraction in a magnetic quadrupole field. We present the resonance energies and decay 
widths for different values of the field gradient. The resulting spectrum consist of several 
well-separated parts. Concerning the resonance positions one can distinguish three regimes, 
each of which reveals individual characteristics: the weak, the intermediate, and the strong 
gradient regimes. In the weak gradient regime, the Zeeman term Hb is very small compared 
with the hyperfme interaction Hhf and only slightly perturbs the zero-field eigenstates of 
H. In this case the atom, being primarily in its hyperfme ground state F = I — |, remains in 
this manifold and behaves approximately like a neutral particle of spin F with the ^-factor 

Q „, F(F + 1)- 1(1+1) +S(S+1) 

9f-9e 2F(F + 1) • [iZ) 

This regime is also called the Zeeman regime. In the intermediate gradient regime the 
Zeeman and the hyperfme interactions are of the same order of magnitude, and the atom 
in the ground state may represent a significant admixture of different hyperfine states F. 
Finally in the strong gradient regime the Zeeman term dominates the hyperfine energies 
at least for sufficiently large distances from the coordinate origin. In this case the spin 
component of the electron along the local direction of the magnetic field is almost conserved 
and as we will see, the resonance positions in the complex energy plane are grouped according 
to different values of its quantum number m s . This regime is also called the Paschen-Back 
Regime. 

Knowing the resonance eigenfunctions of the complex-scaled Hamiltonian one can cal- 



culate corresponding expectation values within the generalized inner c-product [20|. The 
expectation value which is obtained in this way is in general complex. The real part repre- 
sents the average value, whereas the imaginary part can be interpreted as the uncertainty 
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of our observable in a measurement when the system is prepared in the corresponding res- 



onance state 



20| | . We will analyze the average values of the components of the spins which 



point along the local direction of the field as a function of the energy. This enable us to 
explain different sets of resonances in each regime. We also discuss in this section the de- 
pendence of the decay width of a resonance state on its angular momentum as well as on 
the field gradient. 



A. Resonance Positions in the Zeeman Regime 

For large values of (3 ((3 > 200) we are in the Zeeman regime. In Fig. [Tfa-c) we present the 
energies and decay widths for (3 = 1000 and m = 10, for an atom with nuclear spin I — | 
being in its hyperfine ground state. The resonances are localized in the negative energy 
region and their distribution consists of three well-separated parts. They can be classified 
according to the expectation value mp = i?e((F-B/|B|)) which is the projection of the total 
spin onto the local direction of the magnetic field. There are long-lived resonances which 
can be identified with mp = — 1, and two sets of resonances with shorter lifetimes whose 
spin projections are mp = 0, +1. The values of mp are shown in the panel d. For mp = — 1 
one observes the resonances to be located on lines with similar slopes covering the area of a 
right triangle in the e — T plane. For mp = and = +1 one immediately notices that 
both sets are arranged on lines with an infinity and a negative slope respectively. In Figflje) 
we present the decay width of the energetically lowest state of the long-lived resonances 
(the resonances with mp = —1) as a function of the angular momentum m. The decay 
width decreases exponentially with increasing value of m. We also present the decay width 
of the energetically highest state of the two sets of short-lived resonances as a function of 
the quantum number m. Unlike the case mp = —1, the decay width of resonances with 
mp = +1 and mp = increases with increasing angular momentum, (for an explanation 
of this behavior see [7]). FigOJf) shows the expectation value of the squared total spin, 
Re((F 2 )) as a function of the energy for the same parameter values. For all resonance 
states, the value is approximately +2 which corresponds to F — 1 i.e. the atom behaves like 
a spin-1 particle. 

In Figj2] we present the lifetime of the energetically lowest long-lived resonance state 
(mp = —1) for m = 25 as a function of the parameter f3. The lifetime decreases as the 
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FIG. 1: Decay width and energies of the resonances for [3 = 1000, m = 10 for (a) mp = 1, (b) 
rnj? = and (c) = -1. (d) Expectation values of the total spin component along the local 
direction of the magnetic field for the three sets of resonances, (e) (■) The decay width of the 
energetically lowest state with mp = — 1, (A) the decay width of the energetically highest state 
with in p = 0, (*) the decay width of the energetically highest state with mp = +1 as a function 
of the quantum number m. (f) The expectation value of the total spin squared as a function of the 
energy. 



8 



10 2 



10 



■2-10° 



10" 



10" 



10 



10 10 
p [a.u.] 



10" 



FIG. 2: The lifetime of the energetically lowest long-lived resonance state for m = 25 as a function 
of the parameter (5. 



hyperfine parameter (3 decreases, i.e. the gradient field b increases. Performing a line fit we 
find the dependence T[s] » 8.3 x 10~ 4 ^ 2 / 3 [a.u.]. 



B. Resonance Positions in the Intermediate Regime 

Let us now focus on the intermediate regime covering the values 0.2 < (3 < 200, where 
the Zeeman and hyperfine interactions become comparable. We observe three different types 
of behavior. For 40 < (3 < 200, an atom being in the hyperfine ground state, still behaves 
approximately like a spin-1 particle. The resonance spectrum is localized in the negative 
energy domain and consists of three well-separated parts to which we can assign the values 
mp = —1, 0, +1. For the two short-lived states, following the same reasoning as above, the 
width increases when m increases, while for the long-lived states it decreases. However in 
contrast to the Zeeman regime (see Fig [2]), the lifetime of the long-lived states (mp = — 1) 
increases when the hyperfine parameter (3 decreases, i.e. the field gradient b increases (see 
Figj3]). For (3 < 40 this value decreases very rapidly when (3 decreases. In this case the 
resonance states which correspond to higher hyperfine levels F = 2 are more stable, and 
are localized in the positive region of the spectrum. FigJU^a) shows the decay widths and 
energies of the resonances possessing positive energies for /3 = 19. We observe four different 
curved and triangular shaped regions with distinct classes of resonances. In FigJH(b) we 
present the expectation values of the squared total spin, Re((F 2 )), of these resonances as 
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FIG. 3: The lifetime of the energetically lowest long-lived resonance state for m = 25 as a function 
of the parameter (5. 



a function of the energy. The values range from approximately 5.5 to 6.0, correspond to 
F 2.0. FigJH(c) shows the expectation values of the total spin component along the local 
direction of the field for the resonances of FigJU(a). The resonances are divided into four 
well-separated parts corresponding to rrip = +2, +1, 0, —1. Resonances with rrtp = —2 have 
very short lifetime and are not shown. For j3 < 3 the non- conservation of F becomes even 
more explicit, see FigJS^a) for (3=1. In Fig. E(b,c) we present the respective decay widths 
and energies. The resonances which are localized in the negative energy region have short 
lifetimes and are divided into four subgroups. Each group lies on a line with negative slope, 
the four slopes being very similar. Moreover, the resonances investigated form subgroups 
on these lines. The positive energy resonances have much longer lifetimes and cover an 
area of approximately triangular shape in an irregular manner, i.e. no pattern is visible. In 
FigE^d) the corresponding expectation values of the electronic spin component along the 
local direction of the magnetic field, Re({S • B/|B|)) is shown. For the negative energy 
resonances this value is approximately -0.5 indicating that the spin is aligned opposite to 
the local direction of the magnetic field, while for positive energies the spin is parallel to the 
magnetic field. Figj5|e) shows the expectation values of the nuclear spin component along 
the local direction of the field. In the negative energy domain of the spectrum the pattern 
divides into four parts with values mj = —§,—§,+§,+§, while in the positive energy region, 
the pattern is strongly disturbed and mj is not conserved. Figj5]^f) presents the lifetime of 
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FIG. 4: (a) Decay width and energies of the long-lived resonances with higher hyperfine quantum 
number F = 2 for [5 = 19 and m = 25. (b) The expectation value of the total spin squared as a 
function of the energy, (c) Expectation values of the total spin component along the local direction 
of the magnetic field as a function of the energy. 

the energetically lowest long-lived state as a function of (3. This value increases when (5 
decreases, i.e. the gradient field b increases. Our results show that for the short-lived states, 
i.e. the negative energy resonance states, the width increases when m increases, while for 
the long-lived states, i.e. the positive energy resonance states, it decreases. 



C. Resonance Positions in the Hyperfine Paschen-Back Regime 

The hyperfine Paschen-Back regime includes f3 < 0.2. In case the Zeeman term dominates 
the hyperfine energy, it is natural to decompose the Hamiltonian according to H = H + Hi 
where H = |(p 2 + xS x + yS y — 2zS z ) describes the motion of a spin \ particle in the 
magnetic field, and H x = a(xl x + yl y — 2zl z ) + (31 ■ S perturbs the eigenstates of H . The 
spectrum consists of two parts: again we have one set of resonances localized in the negative 
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FIG. 5: (a) The expectation value of the total spin squared as a function of the energy for m = 10 
and P = 1.0. Decay width and energies of the resonances for (b) ms = — |, (c) ms = +|. (d) 
Expectation values of the electronic spin component along the local direction of the magnetic field, 
(e) Expectation values of the nuclear spin component along the local direction of the magnetic 
field, (f) The lifetime of the energetically lowest resonance state of the positive energy domain of 
the spectrum as a function of (3. 
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energy region with short lifetimes and a second set localized in the positive energy domain 
possesses much larger lifetimes. In Figj6|a,b) we present the energies and decay widths 
of resonances for (3 = 0.1 and m = 10. The negative energy resonances are arranged on 
lines with a negative slope. The positive energy resonances cover an area of triangular 
shape, some of them being located on straight lines. Figl^c) shows the corresponding 
expectation value i?e((F 2 )) of the squared total spin. This value is almost 4.5 which again 
indicates that the total spin quantum number F is not conserved. In Figj6]^d) the respective 
expectation value Re((S ■ B/|B|)) of the electronic spin component which points along the 
local direction of the magnetic field is presented. For the positive energy resonances this 
value is approximately +0.5 indicating that the spin is aligned parallel to the local direction 
of the magnetic field while for negative energy resonances it is antiparallel. The nuclear 
spin component along the local field is not conserved, however it is instructive to consider 
its average value, Re((l ■ B/|B|)), as a function of the energy. For large values of m (see 
Figlf^e) for m = 10) the resonances are grouped into four regular subparts for both the 
positive and negative energy domain. 

We remark that the decay widths of the positive energy resonances show again an ex- 
ponential decaying behavior as a function of the quantum number m whereas the decay 
widths of the negative energy resonance states increases with increasing m. The lifetime of 
the energetically lowest long-lived resonances for m = 25 as a function of (3 can be fitted 
according to T[s] « 1.1 x 10 -1 /3 2/3 [a.u.]. 

V. SUMMARY 

We have investigated the resonant quantum properties of an atom with a single ac- 
tive valence electron taking into account its hyperfine structure (J = |) in a 3d magnetic 
quadrupole field. 

We have calculated and analyzed the energies and decay widths of the resonance states 
of the Hamiltonian employing the complex scaling approach and a Sturmian basis set. With 
respect to the resonance position one can distinguish essentially three regimes. In the weak 
gradient regime, the Zeeman term is very small compared with the hyperfine interaction 
and only slightly perturbs the zero-field eigenstates of the Hamiltonian. In this case the 
atom being primarily in its hyperfine ground state behaves approximately like a particle of 
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FIG. 6: Decay width and energies of the resonances for j3 = 0.1, m = 10 for (a) ms = — g, (b) 
ms = +|. (c) The corresponding expectation value of the total spin squared as a function of the 
energy, (d) Expectation values of the electronic spin component along the local direction of the 
magnetic field for the two sets of resonances, (e) Expectation values of the nuclear spin component 
along the local direction of the magnetic field. 

spin 1, and the resonance states are grouped into three well-separated parts, corresponding 
to three different directions of the total spin with respect to the local direction of the field. 
The resonances for which the total spin is antiparallel to the field (m F = —1) possess the 
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largest lifetimes. As the total angular momentum of the atom increases, the decay width of a 
short-lived resonance state increases while for a long-lived resonance state it decreases. Our 
results for the long-lived resonances show that as the scaled hyperfine parameter (5 increases 
so does the lifetime. 

In the intermediate regime the Zeeman and hyperfine interactions are comparable. For a 
sufficiently weak field gradient the atom in the hyperfine ground state behaves approximately 
like a spin-1 particle, and the resonances are arranged somehow similar to the Zeeman regime 
but in contrast to the Zeeman regime its long-lived resonance states become more stable 
when the gradient field increases and/or (3 decreases. For stronger field gradients and/or 
smaller values of the parameter /3 the lifetime decreases rapidly. In this case resonances 
which correspond to the higher hyperfine level, F = 2, are more stable. For even smaller 
values of (3, the non-conservation of F becomes remarkable. In this case the electronic spin 
component along the local direction of the magnetic field is almost conserved. A resonance 
which is localized in the negative energy region of the spectrum possesses a short lifetime. 
In such a state the electronic spin is aligned opposite to the local direction of the magnetic 
field. The positive energy resonances possess a much longer lifetime with their electronic 
spin being parallel to the magnetic field. For a short-lived state the width increases with 
increasing m, while for a long-lived state it decreases. For the latter our results show that 
the lifetime of this state increases with decreasing values of (3. 

In the strong gradient regime, the Zeeman term dominates the hyperfine energy. Here the 
component of the electronic spin along the local direction of the field is almost conserved. 
The energy spectrum is arranged into two disconnected parts each of which contains exclu- 
sively resonance states with negative and positive energies. The former are of short-lived 
character and the electronic spin is antiparallel to the local direction of the magnetic field. 
The latter possess much longer lifetimes and the spin of the electron is parallel to the field. 
As the total angular momentum of the atom increases, the decay width of a negative energy 
resonance state increases while for a positive energy resonance state it decreases. Similar to 
the weak gradient regime the lifetime of a positive energy resonance state, decreases with 
decreasing f3. 
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